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ABSTRACT 

Low-mass stars in the He-core-burning phase (HeCB) play a major role in stellar, 
galactic, and extragalactic astrophysics. The ability to predict accurately the proper¬ 
ties of these stars, however, depends on our understanding of convection, which remains 
one of the key open questions in stellar modelling. We argue that the combination of 
the luminosity of the AGB bump (AGBb) and the period spacing of gravity modes 
(ABi) during the BeGB phase, provides us with a decisive test to discriminate between 
competing models of these stars. We use the MESA, BaSTI, and PARSEC stellar evo¬ 
lution codes to model a typical giant star observed by Kepler. We explore how various 
near-core-mixing scenarios affect the predictions of the above-mentioned constraints, 
and we find that ABi depends strongly on the prescription adopted. Moreover we 
show that the detailed behaviour of ABi shows the signature of sharp variations in the 
Brunt-Vaisala frequency, which could potentially give additional information about 
near-core features. We find evidence for the AGBb among Kepler targets, and a first 
comparison with observations shows that, even if standard models are able to repro¬ 
duce the luminosity distribution, no standard model can account for satisfactorily the 
period spacing of BeCB stars. Our analysis allows us to outline a candidate model to 
describe simultaneously the two observed distributions: a model with a moderate over¬ 
shooting region characterized by an adiabatic thermal stratification. This prescription 
will be tested in the future on cluster stars, to limit possible observational biases. 

Key words: stars: evolution - asteroseismology - stars: low-mass - stars: interiors - 
stars: late-type 


1 INTRODUCTION 

The helium-core-burning (HeCB) phase of low-mass stars 
affects many aspects of their subsequent evolution, and has 
therefore been the subject of numerous studies in the lit- 
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erature. In particular, the red clump (RC) is a well-known 
feature in the Hertzsprung-Russell diagram - as well as the 
color-magnitude diagram - of simple and composite stellar 
populations, and is associated to the low-mass, metal-rich 
stars in the He-core-burning phase (Cannon 1970). The RC 
plays a key role in many helds of astrophysics: its luminos¬ 
ity, for instance, can be used as a distance and age indicator 
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of clusters and nearby galaxies, while the observed chemical 
composition of its members is useful to investigate the chem¬ 
ical evolntion of galaxies (e.g., see Girardi & Salaris 2001; 
Catelan 2009; Nidever et al. 2014). 

Current models of the internal structure and evolution 
of such stars, however, still suffer from systematic uncertain¬ 
ties which are due in most cases to our limited understanding 
of the physical processes in stellar conditions (energy trans¬ 
port and nuclear processes). In particular, predictions of stel¬ 
lar lifetimes in the HeCB phase are strongly dependent on 
the (poorly constrained) amount of mixing applied beyond 
boundaries of convective regions and on the adopted defini¬ 
tion (and correct implementation) of such boundaries (e.g. 
see Chiosi 2007; Castellani et al. 1971a; Gabriel et al. 2014; 
Bressan et al. 2015). Crucially, our ability to test models has 
been limited so far by the lack of observational constraints 
which are specific to the internal structure of evolved stars. 

Asteroseismology of thousands of red giants observed 
by CoRoT (Baglin et al. 2006) and Kepler (Borucki et al. 
2010) has changed the situation. We can now use the pul¬ 
sation frequencies to place tight constraints not only on the 
fundamental stellar properties, but also to probe their inter¬ 
nal structure (see e.g. Chaplin & Miglio 2013 and references 
therein). In particular, as presented in Montalban et al. 
(2013), the frequencies of oscillation modes detected in 
HeCB stars are sensitive diagnostics of the chemical and 
thermal stratification of the energy-generating core, provid¬ 
ing us with a novel and independent constraint, which is 
specific to the core structure of HeCB stars. In this paper 
we argue that seismic constraints, when used in conjunc¬ 
tion with well-known classical indicators (in particular the 
luminosity of the AGB bump, see below), provide us with a 
decisive test to discriminate competing models of different 
near-core mixing schemes adopted during the He-burning 
phase. We also show how current standard models (see Sec. 
3) from different stellar codes do not match these constraints 
simult aneously. 

The paper is organised as follows: we start by review¬ 
ing currently available indicators of core mixing in He-core- 
burning stars (Sect. 2). We describe in Sec. 3 the set of stel¬ 
lar models we use, while in Sec. 4 we explore the impact of 
several commonly adopted prescriptions for near-core mix¬ 
ing on the observables, with emphasis on asteroseismic con¬ 
straints (Sec. 5). We then compare in Sec. 6 our predictions 
with asteroseismic constraints as published in recent papers 
(Mosser et al. 2014; Pinsonneault et al. 2014). Finally, we 
present a summary and future prospects in Sec. 7. 


2 OBSERVATIONAL CONSTRAINTS ON 
HE-CORE-BURNING MODELS 

One of the main observables used to constrain the mass of 
the fully mixed core during the HeCB phase is the R 2 ratio 
(Buonanno et al. 1985). R 2 is defined as the ratio between 
the number of early Asymptotic Giant Branch (eAGB) 
Horizontal Branch (HB) stars in simple stellar populations 
(chemically homogeneous and coeval stars) and is directly 
connected to the lifetime of the two phases { R2 'Tagb/'Tib)- 
The value of R 2 is affected by the core mixing during the 
HeCB phase, as discussed in, e.g., Bressan et al. (1986), 
Caputo et al. (1989). The use of R 2 has hitherto been lim¬ 


ited to stars in clusters, due to the small number of held 
stars with accurate distance available, and due to the inher¬ 
ent complication of dealing with an ensemble of stars with a 
spread in age and chemical composition (which also hampers 
robust inferences on their evolutionary state). 

Another important observable related to the HeCB and 
AGB evolution of low- and intermediate-mass stars is the 
luminosity of the AGB bump (AGBb, see e.g. the review by 
Catelan 2007). Similarly to the bump in RGB, the AGBb 
is a feature of the eAGB that consists in three passages 
of the evolutionary track in a small interval of luminos¬ 
ity and causes a local increment of stars in the luminosity 
distribution of a stellar population. The occurrence of the 
AGB bump is connected to the formation of the helium¬ 
burning shell (see e.g. Cassisi & Salaris 2013, paragraph 
5.2). The luminosity at which it occurs depends on the lo¬ 
cation in mass of the He shell at its ignition, hence it is 
determined by the maximum extension of the mixed core 
during the HeCB phase. While the AGB bump had been 
highlighted in stellar evolutionary tracks a few decades ago 
(see, for instance Caputo et al. 1978), its hrst identihca- 
tion as a distinct observational feature in galaxies was re¬ 
ported by Gallart (1998). The AGBb is nowadays observed 
in a large number clusters and nearby galaxies (e.g., see 
Alcock et al. 1997; Ferraro et al. 1999; Beccari et al. 2006; 
Dalcanton et al. 2012). A further important property of the 
AGBb is the weak dependence of the luminosity ratio be¬ 
tween the RC and the AGBb on the metallicity and initial 
helium abundance (Castellani et al. 1991; Bono et al. 1995). 

Asteroseismology of white dwarfs may also provide us 
with observational constraints to test models in the HeCB 
and AGB phase. The C/O profile of a white dwarf at the 
beginning of the cooling sequence corresponds to that of the 
stellar core at the end of the AGB, hence it is largely de¬ 
termined by the combined action of chemical mixing and 
nuclear burning during the HeCB phase. Using pulsation 
modes detected in WDs, Metcalfe et al. (2002) found a dis¬ 
crepancy in the central oxygen-to-carbon ratio between stel¬ 
lar models and the value inferred from seismic data. They as¬ 
cribed this difference to an underestimation of *^C( 0 !, 7 )^®O 
cross section, however, as pointed out by Straniero et al. 
(2003), the final C/O ratio in models also depends on the 
amount of mixing applied in the convective regions beyond 
the Schwarzschild border, and on the adopted dehnition of 
such boundaries. 

While providing valuable information, currently avail¬ 
able tests cannot be used to discriminate competing models. 
In this paper we propose a way forward, which is based on 
the combination of said constraints and the more direct di¬ 
agnostics of the conditions in the core provided by non-radial 
modes observed in HeCB stars. 


3 STELLAR MODELS 

In our exploratory analysis we consider models with 
M=1.5Mq and solar chemical composition. This set of 
parameters are both typical of the sample of RC gi¬ 
ants observed by Kepler (see e.g. Pinsonneault et al. 2014; 
Mosser et al. 2014), and similar to those of giants in the 
cluster NGC6819 (see Basu et al. 2011; Miglio et al. 2012; 
Sandquist et al. 2013). Models in HeCB phase are computed 
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using three different stellar structure and evolution codes 
(MESA, BaSTI, and PARSEC) and several assumptions 
about near-core mixing. 

A first set of models presented in this study is computed 
using the Modules for Experiments in Stellar Astrophysics 
code (MESA, Paxton et al. 2011). The choice of parameters 
and of the relevant physics is the following: 

• the initial mass is 1.5Mq, and no mass loss is considered 
during the RGB; 

• the mass fraction of heavy elements (Z = 0.01756), the 
initial helium abundance {Y = 0.26556), and the mixing- 
length parameter (o!„it = 1.69278) are those resulting from 
a solar-calibrated model without diffusion; 

• the partition of heavy elements is taken from 
Grevesse & Noels (1993); 

• the nuclear reaction rates are taken from Angulo et al. 
(1999); 

• during the He-burning stage one has to properly ac¬ 
count for the C and O abundance changes in the stellar 
matter and their effect on radiative opacities, hence we use 
opacity tables accounting for various levels of C and O en¬ 
hancements (Type 2 opacity tables, Iglesias & Rogers 1996); 

• the evolution includes the pre-main sequence phase and 
is computed until the first thermal pulse along the AGB 
(TP-AGB); 

• in models computed with overshooting, the overshoot¬ 
ing scheme adopted is a step function in which the parameter 
O^ovHe indicates the extension of the mixed region from the 
Schwarzschild border in units of the local pressure scale- 
height. The overshooting region is instantaneously mixed 
(Maeder 1975); 

• we consider two scenarios for the temperature gradi¬ 
ent (Vj-) in the overshooting region, according to the defini¬ 
tion given in Zahn (1991); V 7 -=Vi- (overshooting) or V 7 -=Va(j 
(penetrative convection), where and follow the usual 
notation for the radiative and adiabatic temperature gradi¬ 
ent; 

• during the main-sequence phase we adopt an overshoot¬ 
ing parameter aovH=0.2 (Nevertheless, the choice of ttovn 
does not impact on the structure of a low-mass HeGB star). 

We compute the following MESA models with different 
schemes for convective mixing: 

BS The “bare Schwarzschild” model (BS, to follow the 
notation in Straniero et al. 2003). In the BS model, the 
Schwarzschild criterion is applied on the radiative side of 
the convective boundary which implies that the convective 
core is not allowed to grow, leading to a discontinuity in the 
chemical composition (hence in the radiative gradient) at 
the border of the convective core (Figure 1 panel a). 

A recent paper by Gabriel et al. (2014) showed how this 
implementation of the Schwarzschild criterion leads to an in¬ 
consistent location of the convective border, as the convec¬ 
tive luminosity is non-zero there (Vr > Vad at the inner side of 
the convective boundary). The book of Schwarzschild (1958) 
and the analysis presented in Gastellani et al. (1971a) also 
led to similar conclusions. The latter consider the BS convec¬ 
tive border as in an unstable equilibrium, in the sense that 
by extending outwards the convective core by an arbitrar¬ 
ily small quantity, Vr at the new border may be larger than 
Vad, hence, the region is convectively unstable according to 


the Schwarzschild criterion. To find a stable border, the con¬ 
vective core must be extended until radiative and adiabatic 
gradient become equal (induced overshooting). However, as 
the HeGB proceeds (Y^ < 0.69) this scheme generates a local 
minimum in Vj- within the mixed core. The addition of more 
radiative layers surrounding the increasingly larger convec¬ 
tive core will decrease Vj to the value of Vad fhe location 
of the minimum. This creates a separate convective region 
in layers beyond the location of this minimum, where V^ is 
still larger than Vad. The treatment of this external convec¬ 
tive region is problematic. Full mixing between layers inside 
the minimum of V^ and the external convective shell can¬ 
not happen, because otherwise the minimum of Vj would 
decrease below the local value of Vad with the contradiction 
of a fully mixed region where V,- is however locally lower 
than Vad. A solution to this problem is the formation of a 
partially mixed -semiconvective- region between the min¬ 
imum of Vr and the outer radiative zone (Gastellani et al. 
1971b). This is usually treated with dedicated algorithms 
that allow for partial chemical mixing to satisfy Vj- = Vad in 
this region, with the consequence of creating a smooth gra¬ 
dient of chemical composition before the sharp discontinuity 
due to the HeGB. Although we discourage to adopt the bare 
Schwarzschild model, we use this mixing scheme in order to 
compare our models with those described in the literature 
(e.g. Straniero et al. 2003; Kippenhahn et al. 2013). 

HOV High overshooting model with CtovHe = 1 (HOV). 

HPC Penetrative Gonvection model with a high over¬ 
shooting parameter ttovHe = 1 (HPG). 

For low and moderate values of ttcvHe, overshooting 
models (as well as penetrative convection) might develop 
a semiconvective region similar to the case presented by 
Gastellani et al. (1971b). For high values of the overshoot¬ 
ing parameter (e.g HOV and HPG) the extra mixed region 
becomes large enough to prevent the formation of a semicon¬ 
vective region (see e.g. Bressan et al. 1986; Straniero et al. 
2003). Moreover, in these cases the determination of the fully 
mixed region requires the application of the Schwarzschild 
criterion in layers where there is no chemical composition 
gradient/discontinuity, which greatly simplifies the numeri¬ 
cal implementation of the convective-instability criterion. 

In order to extend our study of the mixing-schemes com¬ 
monly adopted, we also consider stellar models from other 
evolution codes (BaSTI and PARSEG). To compute those 
models we used, where possible, the same physical inputs 
adopted in MESA. 

BaSTI-SC We have also computed models with the 
BaSTI code (A Bag of Stellar Tracks and Isochrones, 
Pietrinferni et al. 2004, 2006, 2013). The BaSTI model has 
M = 1.5M0, computed for Z = 0.0176, Y = 0.266, and otnit = 
1.69 with Grevesse & Noels (1993) heavy element partition. 
The input physics relevant to this analysis is the same as in 
the MESA calculations but for the *^C(0!,y)*®O reaction rate 
that comes from Kunz et al. (2002). Gore mixing during the 
HeGB stage, induced overshooting and semiconvection have 
been taken into account by adopting the numerical scheme 
firstly introduced by Gastellani et al. (1985) and previously 
described. Indeed, a semiconvective region starts to appear 
when the central abundance by mass of He is 0.7. The oc¬ 
currence of breathing pulses - which appear when the central 
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abundance of He drops below of 0.12 - is inhibited by impos¬ 
ing that the abundance of Helium in the core is not allowed 
to increase at each time step. The evolution starts from the 
pre-main sequence to the RGB-tip, and is resumed at the 
start of the HeCB after the core electron degeneracy has 
been removed (the helium-flash evolution is not computed), 
and 3% of carbon has been produced during the He-flash. 
No core convective overshoot during the central H-burning 
stage and no mass loss during the RGB phase are taken in 
account. 

PARSEC-LOV We also computed models using 
PARSEC (PAdova & TRieste Stellar Evolution Gode, 
Bressan et al. 2012, 2013) for M = I.SMq, Z = 0.017, and 
Y = 0.279. We briefly summarize below the main input 
physics adopted for the current models. The nuclear reac¬ 
tion rates and corresponding 2-values are the recommended 
values in the JINA reaclib database (Cyburt et al. 2010). 
The high-temperature opacities, 4.2 < log(r /K) < 8.7, are 
provided by the Opacity Project At Livermore (OPAL) 
team (Iglesias & Rogers 1996, and references therein) and 
the low-temperature opacities, 3.2 < log(r /K) <4.1, are from 
jESOPUS^ tool (Marigo & Aringer 2009). The equation of 
state is computed with the FreeEOS code (A.W. Irwin^). 
The heavy element partition is from Caffau et al. (2011). 
The MLT parameter, = 1.74, is calibrated on the 
solar model accounting for element diffusion. The standard 
mixing scheme in PARSEG considers an overshooting 
parameter of A; = 0.5 across the formal Schwarzschild 
border, which means about 0.25Hp above it, and a radiative 
thermal stratification of the extra-mixed region. Similar to 
BaSTI, the evolution starts from the pre-main sequence, 
stops at the helium flash and restarts after the core electron 
degeneracy has been removed after the flash. The amount 
of carbon consumed to remove the electron degeneracy is 
computed from the variation of the gravitational binding 
energy of the core during the flash. During the HeCB phase, 
besides accounting for core overshooting, the code may 
deal with residual semi-convective instabilities using the 
Schwarzschild criterion, and suppresses possible breathing 
pulses of convection. 

PARSEC-MPC We use PARSEC to compute a modi¬ 
fied track in which we extend the mixed He core by pene¬ 
trative convection (adiabatic stratification) with a moderate 
value of overshooting (A^ = 1.0, i.e. ttovue ~ 0.5). 

ES Finally, we include in the analysis a model obtained 
with MESA in which we include an external routine to im¬ 
plement a similar prescription as in BaSTI-SC. This routine 
consists in the following steps: 

1. at the beginning of each timestep we set the posi¬ 
tion of the convective-core boundary and fully mix the 
convective region according to the Schwarzschild criterion 
as implemented in MESA. 

2. we let the code calculate the burning during the 
timestep. 

3. at the end of the timestep, we check whether the 
radiative gradient (with the new composition determined 
by the burning) at the convective border is higher than 
the adiabatic gradient. 


' http://stev.oapd.inaf.it/aesopus 
^ http://freeeos.sourceforge.net/ 



Figure 1. Ratio between the radiative and adiabatic temperature 
gradients as a function of mass in MESA HeCB models computed 
using the BS (panel a) and ES (panel b) scheme. The latter leads 
to a self-consistent definition of the convective boundary (follow¬ 
ing the Schwarzschild criterion), as the adiabatic and radiative 
gradients are equal at the boundary of the core. Prom bottom to 
top the lines are for a sequence of models at the stages of L, = 0-93, 
0.87, 0.81, 0.74, and 0.70. 

4a. if this is the case, we restart the timestep from point 
1, but extending the core boundary by one mesh. 

4b. if not, the equality of the gradients has been 
achieved and the code can continue to the next timestep 
(point 1). 

As a result of this algorithm the convective core extends 
naturally (Fig. 1 panel b). We call this model Extrapolated 
Schwarzschild (ES). The internal structure and seismic pre¬ 
dictions of ES (see Sec. 5.1 and Figures 3 and 4) are remark¬ 
ably similar to the BaSTI-SC. Differently from BaSTI and 
PARSEC, however, a routine to treat semiconvection is still 
to be developed in MESA, hence we stop the evolution of 
the ES model at Yc ~ 0.69. 


4 PREDICTED STELLAR PROPERTIES OF 

MODELS IN THE HECB AND AGB PHASE 

We now compare the properties of the series of models pre¬ 
sented in Section 3, with particular emphasis on those which 
can be tested via a direct comparisons to observations. We 
start by comparing predictions of non-seismic observables 
(see Sec. 2), while in Sec. 5 we focus on seismic diagnostics. 

The main consequence of applying different mixing 
schemes in the HeCB phase is to vary the core mass under¬ 
going convective mixing. We refer to convective-core mass 
(Mcc) as the core mass in which Vj > (formal 

Schwarzschild core), while the mixed-core mass (A/mc) in¬ 
cludes also the extra-mixing region (e.g convective core -I- 
overshooting and/or semiconvection). The mixed-core mass 
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Figure 2. Convective-core (upper panel), mixed-core mass (mid¬ 
dle panel) and total luminosity (lower panel) as a function of 
time for the BS, HOV, HPC, BaSTI-SC, PARSEC-LOV, and 
PARSEC-MPC models from the start of the HeCB phase up to 
the first AGB-TP. Both HPC and HOV models show micropulses. 

also indicates the location in the stellar structure of the dis¬ 
continuity in the chemical composition due to HeCB. In the 
ES models (in the domain we were able to explore) and BS 
models, corresponds to Mmc since no extra mixing is 
introduced. This is also true for the BaSTI-SC model when 
Kc > 0.7, i.e. where semiconvection has not appeared yet, and 
during the last stages of the HeCB, where the convective core 
grows rapidly and the size of the semiconvective region is re¬ 
duced (see Figure 2, upper panel). In HOV, HPC, PARSEC- 
LOV, and PARSEC-MPC, and Mmc are distinct from the 
beginning to the end of the HeCB phase (Figure 2, upper 
and middle panels). Although penetrative convection models 
are expected to have smaller Mcc and M^c than overshooting 
models (see Godart 2007; Noels et al. 2010, in the case of 
massive main-sequence stars), max(Mmc) in HOV and HPC 
is very similar, providing a similar LAGBb (see Figure 2, up¬ 
per and lower panels). 

4.1 Luminosity and duration of the HeCB and 
eAGB phases 

An obvious effect of increasing Mmc is to increase the dura¬ 
tion of the HeCB phase (thbCb)- In models with larger Mmc 
more fuel is available for the triple-a and *^C( 0 !, 7 )*®O nu¬ 
clear reactions, increasing the time necessary to deplete all 
the helium in the core, i.e. the lifetime of the HeCB phase. 
In HOV and HPC models, the HeCB lifetime increases by 
about 40% compared to the BS model, while in BaSTI-SC 
and PARSEC-LOV by about 20%. However, if we look at 
the lifespans of the single phases, the increase in duration 


of the HeCB phase corresponds to a (non-linearly propor¬ 
tional) decrease of the duration of the AGB phase, since 
the formation of the He shell takes place closer to the H 
shell. The BS model has shorter lifetime with respect to 
the HPC and HOV models, however, it has a longer AGB 
phase. The BaSTI-SC, PARSEC-LOV, and PARSEC-MPC 
tracks have a behaviour which is in between the BS and 
the HPC/HOV models. In general we expect that models 
with increasing Mmc produce a longer HeCB phase and a 
less populated early-AGB, combined with a decreasing R 2 
factor from BS models to the HOV and HPC models. These 
results are reported in Table 1. 

The stellar luminosity is also affected by the mixing 
scheme adopted during HeCB (Figure 2, lower panel), and 
LAGBb increases when extra mixing is added, while it remains 
almost unchanged when comparing HOV and HPC tracks 
(Table 1). The maximum extension of (which is very 
similar in HOV and HPC models) corresponds to the inner 
border of the He-shell at its ignition and it determines the 
^AGBb- 

Models with larger cores (HOV and HPC) show mi¬ 
cropulses (Mazzitelli & D’Antona 1986), i.e. a series of sec¬ 
ular instabilities that may occur during the formation of the 
He-burning shell. Micropulses appear after the maximum in 
luminosity of the AGB bump if the nuclear reactions in the 
core stop before the shell is ignited in “thin” conditions (see 
Schwarzschild & Harm 1965). To support this statement we 
have looked at the central helium abundance and the nu¬ 
clear energy generation rate at the start of the AGB bump. 
The values found indicate that in the HPG and HOV mod¬ 
els the contribution of HeCB to the luminosity is negligible, 
since helium is almost completely depleted. On the other 
hand, in the BS model core-He burning still contributes sig¬ 
nificantly to the luminosity. We note that for micropulses 
to be resolved by a stellar evolution code, a small timestep 
is needed in the numerical simulations: each pulse lasts for 
about 5 • 10^ yr and the duration of the entire phenomenon 
is about 2 million years. 

The study of micropulses goes beyond the purpose of 
this paper and for detailed explanations we refer the reader 
to Mazzitelli & D’Antona (1986) and the more recent work 
by Gautschy & Althaus (2007). Also notable of citation are 
the papers by Bono et al. (1997) and Sweigart et al. (2000) 
where the connection between “gravonuclear instabilities” 
and micropulses is investigated. 


5 ASTEROSEISMIC SIGNATURES 

5.1 Asymptotic gravity-mode period spacing 

We have shown that there are no significant differences be¬ 
tween models with similar Mmc when considering the lu¬ 
minosity as a function of time (see e.g. HOV and HPC in 
Fig. 2). These models, however, have very distinct seismic 
properties. Montalban et al. (2013) showed that extending 
the adiabatically stratified central region leads to a larger 
value of the gravity-mode period spacing. The reason why 
the asymptotic period spacing of gravity modes (AIIi) dur¬ 
ing the HeCB phase is a sensitive probe of the temperature 
stratification of near-core regions is directly related to the 
behaviour of the Brunt-Vaisala frequency N, and its relation 
with AHi. 
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Table 1. Non-seismic properties of I.S-Mq models computed adopting different mixing schemes. 



MESA 

BS 

MESA 

HOV 

MESA 

HPC 

MESA 

ES 

BaSTI 

SC 

PARSEC 

LOV 

PARSEC 

MPC 

timespan HeCB phase (thcCb) [Myr] 

60.4 

126.2 

131.1 

— 

100.2 

100.6 

107.0 

R 2 ratio (|ja) 

0.71 

0.06 

0.06 

— 

0.12 

0.15 

0.12 

Mcc at Yc ~ 0.7 [Mq] 

0.100 

0.113 

0.084 

0.162 

0.165 

0.140 

0.107 

Mpic at Yc ^ 0.7 [Mq] 

0.100 

0.292 

0.242 

0.162 

0.165 

0.187 

0.206 

maximum extension Mcc [M©] 

0.100 

0.130 

0.114 

— 

0.243 

0.199 

0.158 

maximum extension Mmc [M©] 

0.100 

0.305 

0.302 

— 

0.243 

0.251 

0.262 


In the stellar interior, N depends on the local tempera¬ 
ture and chemical composition gradients: 

(I) 

where = dln/t/dlnP, Xt = ('?ln/’/dlnr)p^^, xp = 
{dlnP/d\np)Tn, and Xp = {d\nP/d\np,)Tp. In a fully 
mixed region, Vp is null, therefore the profile is directly 
proportional to the difference between V 7 and Vaj (Figure 
4a). 

In the extra-mixed region of models with overshooting 
(e.g. HOV and PARSEC-LOV) Vr < V^d, therefore 

oc Vad — Vj- = Vad — Vr > 0, 

while in the corresponding region of a penetrative convection 
model (e.g. HPC and PARSEC-MPC) V 7 = Vad, hence 

‘Xi Vad ~ V7 =0. 


The asymptotic period spacing of gravity modes is re¬ 
lated to the Brunt-Vaisala frequency according to the rela¬ 
tion (Tassoul 1980): 


AH; 


\/l(l + 1) \J''i r J 


( 2 ) 


where ^2 and rj are the boundaries of the g-mode cav¬ 
ity (i.e. where or = A^). Consequently, HeCB penetrative- 
convection models have higher values of AIli compared to 
overshooting models (Figure 3). 

A similar effect can be found if we compare models 
with increased mixed-core size. Bare-Schwarzschild models, 
in fact, have lower AIIi compared all other models, followed 
by models with induced overshooting (BaSTI-SC and ES), 
then high overshooting model, and finally high penetrative 
convection models (see Fig. 4a). 

Independently from the convective-mixing scheme 
adopted (BS, HOV and HPC), the period spacing of mod¬ 
els in the HeCB is higher compared to that on the RGB at 
the same luminosity (AHiRgg ~ 60 —50 s), while after the 
early-AGB phase AHi decreases to similar or smaller values 
(Montalban & Noels 2013). 

It is interesting to notice that if models present multiple 
gravity-mode cavities, we do not expect them to show a reg¬ 
ular period spacing. This is the case of HOV and HPC mod¬ 
els in the in the post-HeCB phase during the micropulses. 
In fact, the He shell can experience convection within each 



Figure 3. Period spacing against the central helium mass frac¬ 
tion during the HeCB phase for the model computed. Bare- 
Schwarzschild models present the smallest period spacing, while 
the penetrative convection models have the largest values. More¬ 
over penetrative convective models start with larger period spac¬ 
ing. 


pulse, in conjunction with the maximum of the nuclear en¬ 
ergy generation. Stars presenting this scenario (if any at all 
exist) may be missed by analysis based upon looking for a 
simple pattern in AIIi (e.g. Mosser et al. 2012). 


5.2 Period spacing of numerically computed 
adiabatic frequencies 

From the observed frequency spectra we can estimate 
AHi based on the detection of modes that have rela¬ 
tively small inertias, i.e. gravity modes that have a signifi¬ 
cant coupling with the low-inertia acoustic modes (see e.g. 
Christensen-Dalsgaard 2012 and reference therein). The fre¬ 
quencies of these mixed gravito-acoustic modes are expected 
from theory to follow a relatively simple pattern (see e.g. 
Unno et al. 1989; Mosser et al. 2012), which can be fitted 
to the observation to estimate AIIi, provided that a suffi¬ 
cient number of modes are detected, and that the analytical 
approximation for the expected pattern of mixed modes is 
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Figure 4. (a): Brunt-Vaisala frequency in the stellar interior of models with ~ 0.7 and different convective-mixing scheme. Figures 
(b-h): Oscillation modes properties for the different convective schemes. Upper panel: mode inertia as a function of the frequency for 
modes with angular degree ^ = 0,1,2. Lower panel: period spacing of numerically computed dipolar-mode frequencies (dots) compared 
with the asymptotic value (solid orange line). 
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Figure 5. Brunt-Vaisala frequency as a function of the nor¬ 
malised buoyancy radius n,/nj( (see text) in a HOV model with 
Yc ^ 0.7. In the rectangle, £ = I period spacing as a function of 
period for the same model. The period spacing of high-order g 
modes (e.g. in the period range 4.7 —5.3 10"^ s) is well described 
by the superposition of the asymptotic AHi (orange line) and a 
component with periodicity An ~ 6. This periodicity indicates (see 
Eq. 3) a sharp-structure variation located at n^/rir ~ 1/A« ~ 0.17, 
which corresponds well with the position of glitch in N. 


accurate (see e.g. Beck et al. 2011; Montalban et al. 2013; 
Mosser et al. 2012; Jiang & Christensen-Dalsgaard 2014). 

It is thus crucial not only to make predictions of RC 
AIli using the approximated expression in Eq. 2, but also 
to compute the spectrum of individual modes, which may 
show interesting departures from the expected approximated 
relation/pattern, and that eventually can be compared with 
the detailed observed frequency spectrum. 

In Figures 4b— f we present the properties of 
adiabatic pulsation modes computed with GYRE 
(Townsend & Teitler 2013), and compare frequency 
spectra of models with the same central helium abundance 
(Tc — 0.7), but computed with different convective-mixing 
schemes. In the upper panel of each figure it is possible to 
see how the inertia (E) of each mode varies in a frequency 
interval around the frequency of maximum oscillations 
power (Vtnax)- Figures 4b— f also show that the characteristic 
asymptotic behaviour of the modes (the constant frequency 
separation for the low-inertia, pressure-dominated modes 
and the constant period spacing for the high-inertia, 
gravity-dominated modes) is a good representation of 
the detailed, numerically computed frequency spectrum. 
Moreover, the asymptotic value of the period spacing clearly 
reflects the differences in the Brunt-Vaisala frequency near 
the core (Fig. 4a and Eq. 2), with the BS model having the 
lowest AIIi, and the HPC model the highest. 


5.3 Signatures of sharp-structure variations in the 
period spacing 

As evinced from Fig. 4a and, more clearly, from Fig. 5, N may 
have sharp changes due e.g. to chemical composition gradi¬ 
ents and/or in the temperature gradient in radiative regions. 
Whether such glitches have a significant impact of the pe¬ 
riod spacing depends on their location, their sharpness, and 


to the typical local wavelength of the gravity modes of in¬ 
terest. 

As described in the literature (see e.g. Brassard et al. 
1992; Miglio et al. 2008; Berthomieu & Provost 1988), the 
signature of a sharp feature in the Brunt-Vaisala frequency 
is a periodic component in the periods of oscillations, and 
therefore in the period spacing, with a periodicity in terms 
of the radial order n given by: 


where the total buoyancy radius is defined as: 

n«' = / (4) 

•/ro ^ 

and local buoyancy radius is 

ny' = / -dr ', (5) 

Jro r 

with ro and R being the inner and outer boundary of the 
g-mode propagation region. The periodicity of the compo¬ 
nents in terms of radial order is therefore an indicator of the 
location of the glitch, expressed in terms of its normalised 
buoyancy radius (fls/TIgiitch)- 

In the models we are focusing on, we notice two main 
glitches. A prominent, yet smooth, glitch due to the FI- 
burning shell and associated with (see e.g. the glitch 
located at nfi/IIr cs 0.7 in Fig. 5). In models at the begin¬ 
ning of the HeCB phase, in which the H-burning shell is 
still very thin, this glitch is sharper and may give rise to 
significant departures from a smooth g-mode period spacing 
(see Bildsten et al. in preparation). In models during most 
of the HeCB phase, however, this glitch does not appear to 
give rise to significant deviations from the asymptotic AIIi 
expected for high-order g modes. 

More interestingly, in the HOV model (see Fig. 4c and 
5) we notice a sharp variation in N which can be well de¬ 
scribed by a step function^. Given the location of this glitch 
Hr/II;- ~ 0.17 (see Fig. 5) we expect a periodic deviation 
from the asymptotic AHi with a periodicity of An ~ 6, where 
n is the radial order of gravity modes (see e.g. Miglio et al. 
2008). This corresponds well (at least in the regions of pure 
g modes, e.g. in the range 4.7 —5.3 10^ s) to the periodicity 
of the component (see inset of Fig. 5). 

While the full description of these glitches is beyond the 
scope of the present paper, we note that departures from the 
simple description of AHi expected from the interaction be¬ 
tween high-order g modes and an acoustic mode (Unno et al. 
1989; Mosser et al. 2012; Jiang & Christensen-Dalsgaard 
2014, e.g., see) provide additional, potentially very sensitive, 
probes of sharp-structure variations in near-core regions dur¬ 
ing the HeCB phase. 


^ At the boundary of the fully mixed radiative (overshoot) re¬ 
gion, the discontinuity in opacity, due to the difference between 
carbon rich mixed layers and He rich surrounding layers, leads to 
a discontinuity in Vj- and thus in N. 
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Figure 6. Hertzsprung-Russell diagram and luminosity distribu¬ 
tion of the APOKASC catalogue (Pinsonneault et al. 2014). The 
sample is selected in a range of mass M= 1 . 3 - 1 . 7 M 0 and metal- 
licity [M/H] =-0.4-0.4. The AGB bump can be identified in the 
right the panel at around logL/L 0 ~ 2.2 (black arrow) while the 
major peak correspond to the RC. 


6 FIRST COMPARISON WITH 

OBSERVATIONS: AGB BUMP AND PERIOD 

SPACING 

6.1 The AGB bump in Kepler red giants 

A catalogue of seismic (Av and Vmax) and spectroscopic 
([M/H] and T/ff) constraints for ~ 1600 Kepler giants 
was recently published by the APOKASC collaboration 
(Pinsonneault et al. 2014). 

We estimate stellar masses, and luminosities of these 
stars by using the so-called direct method, i.e. by combin¬ 
ing T/ff with the seismic radii (estimated using Av and Vmax 
through scaling relations). This method is known to lead 
to less precise estimates of M and R than so-called “grid- 
based” approaches, yet it is less dependent on stellar evo¬ 
lutionary tracks (e.g. see Chaplin & Miglio 2013, and ref¬ 
erences therein) and completely independent of bolomet- 
ric corrections. We select stars from the APOKASC cat¬ 
alogue in a range of mass M =1.3-1.7 Mq and metallicity 
[M/H] =-0.4-0.4. The luminosity function of such stars dis¬ 
plays a peak that is spread over about 4 bins and a has max¬ 
imum at logL/L 0 ~ 2 . 2 : we interpret this peak as a strong 
candidate for the AGB bump (Figure 6 ). This statement is 
supported by the fact that 97.5% of the stars have errors 
smaller than the bin size. We also calculated 1000 realisa¬ 
tions of the observed sample, assuming gaussian errors on 
Av, Vjjiax, Uff (taken from the APOKASC catalogue) and 
found that the properties of the peak in the luminosity func¬ 
tion are not significantly affected. 

As a word of caution we would like to stress that the 
APOKASC catalogue may, however, be affected by target se¬ 
lection biases (Pinsonneault et al. 2014). Although the max¬ 
imum in the observed luminosity distribution does not ap¬ 
pear to be significantly affected by widening the metallicity 
range, we notice that, if we extend the range of masses down 



Figure 7. Hertzsprung-Russell diagram (left panel) and luminos¬ 
ity density distribution (stacked histogram, right panel) of stars 
in the TRILEGAL simulation of the Kepler field for the same 
masses and metallicities as in Figure 6 . The AGB bump is visible 
in the right panel around logL/L 0 ~ 2.2 (black arrow). 


to 1 M 0 , the position of the peak is lowered by 1 bin (0.07 
dex). 

To check whether the AGBb is a feature we expect to be 
able to detect in a composite stellar population we use the 
TRILEGAL code (Girardi et al. 2012) to simulate the galac¬ 
tic population expected in Kepler field. The stellar models 
used are based on Padova tracks (Bressan et al. 2012), with 
an overshooting parameter of A^ = 0.5 (same mixing scheme 
adopted for PARSEC-LOV model). We apply to the syn¬ 
thetic sample the same selection in mass and metallicity as 
in the observed sample. The synthetic population also dis¬ 
plays a well defined peak with a luminosity compatible with 
that of the candidate AGBb in the APOKASC catalogue 
(Figure 7). 

From the simulations we estimate that about half of the 
red-giant stars in that peak belong to the RGB. This means 
that in the case of the full APOKASC catalogue (considering 
the entire range of mass and metallicity) we expect about 
20 — 40 AGBb stars. 

In order to make a first comparison between our models 
and the observations, we generate histograms of the luminos¬ 
ity based on each of our tracks. Results of this comparison 
are shown in Figure 8 . For the sake of clarity in the figure 
we have omitted RGB models. All the models are able to re¬ 
produce reasonably well the position of the RC, taking into 
account also the fact that our models are representative of 
a single-mass, single-metallicity population only, and that 
we have not added the effect of observational uncertainties 
when building the synthetic luminosity function. 

The luminosity of the AGB-bump predicted by BaSTI- 
SC, and PARSEC-MPC are remarkably similar, and in good 
agreement with the candidate AGB-bump luminosity as de¬ 
tected in the observations. While PARSEC-LOV is another 
acceptable model, interestingly, the AGBb predicted by the 
BS, HOV, and HPC models is in clear disagreement with 
the observations, being either too faint (BS) or too bright 
(HOV, HPC). 
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log L/Lq 


Figure 8. Comparison between the observed luminosity density distribution of stars in APOKASC public catalogue selected with the 
same criteria as in Figure 6, and the predictions from the models described in the paper. In the models, RGB stars are omitted. 



AHi [s] 


Figure 9. Comparison between the observed period-spacing distribution of stars classified by Mosser et al. (2014) as RC and AGB and 
the predictions from the models described in the paper. We have selected stars with mass 1.3 <M/M(^ < 1.7. 
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6.2 Period spacing of Kepler RC stars 

Mosser et al. (2014) recently published a catalogue that con¬ 
tains a large sample of red giants in the Kepler Held for which 
g-mode asymptotic period spacings were inferred by mod¬ 
elling the interaction between pressure and gravity modes 
(see Unno et al. 1989; Mosser et al. 2012). Thanks to this 
catalogue we are able to compare the observed period spac¬ 
ing with theoretical predictions. We select RC stars (i.e. stars 
with AITi > 175 s) in a range of mass between 1.3 Mq and 
1.7 Mq, and compare the distribution of observed AIli with 
predictions from our models (see Fig. 9). 

None of the models considered seems to describe sat¬ 
isfactorily the entire observed distribution. The BS model 
can be ruled out, since maxAITi is 80 — 90 s smaller than the 
maximum in the observed data. BASTI-SC, which was one 
of our best candidate models based on the comparison with 
the AGBb luminosity, shows a main peak in the AIli dis¬ 
tribution about 20 — 30s below the observed one and cannot 
describe period spacings higher than 300s . PARSEC-LOV 
has a very similar behaviour to BaSTI-SC and has a max¬ 
imum AIli cs 285 s. The HPC model is able to reproduce a 
large domain of the observed distribution, although its AIli 
appears shifted to higher values with respect to the obser¬ 
vations. Finally, although PARSEC-MPC and HOV cannot 
reproduce the highest observed AITi, their main peaks agrees 
with the observed one within 1-2 bins resolution. 

In the observed distribution of AIli (see Fig. 9) we no¬ 
tice, in addition to the main peak, a secondary peak at 
AIIi ~ 253 s, which still belongs to the RC population. In¬ 
terestingly all our models, but BS, show the presence of this 
second structure. The reason behind it can be deduced sim¬ 
ply following the theoretical tracks in the Kc — AIli diagram 
(see Fig. 3). The tracks cross 3 times a narrow interval of 
period spacing: the first two times at the very beginning of 
the HeCB phase, i.e where AIIi decreases until a local mini¬ 
mum and then starts to rise, while the third passage happens 
during the rapidly decrease that follows the absolute max¬ 
imum (at the end of the HeCB). Since time and Pc are, in 
first approximation, linearly related during the HeCB, the 
third passage is very quick compared to the first and the 
second, therefore the peak is almost entirely populated by 
stars at the beginning of the phase. In the BS model the 
main and the second peak are not resolved since AHi stays 
nearly constant for most of the HeCB phase. 

At this stage we are however limited to a qualitative 
comparison between observed and theoretical distributions 
of AIIi. Such limitation arises from the fact that observa¬ 
tional biases at the target selection stage, and in the de¬ 
termination of the AHi from the power spectra, have not 
been fully explored yet. Moreover, although we have cho¬ 
sen models with a mass and metallicity representative of 
the stars observed by Kepler {[M/H] = —0.07±0.24), and se¬ 
lected stars according to the mass, our synthetic population 
is rather simplistic. As a first test of the effect of changing 
the mass and metallicity, we consider HOV models with dif¬ 
ferent mass (M = 1 Mq) or metallicity (Z = 0.007). We find 
that while the AHi distribution depends little on the mass, 
the effect of reducing significantly the metallicity is to ex¬ 
tend to the range of AIIi by ~ 20 s. This effect needs to be 
taken into account when making quantitative comparisons 


between the observed and theoretically predicted distribu¬ 
tions. 

While these limitations can be partly mitigated by a 
more realistic synthetic population, and by a thorough ex¬ 
amination of selection biases, it is likely that the most robust 
inference will be possible when applying our test to simple 
stellar populations, i.e. to red giants in the old-open clusters 
NGC6791 and NGG6819. 


7 SUMMARY AND FUTURE PROSPECTS 

The ability to predict accurately the properties of He-core- 
burning stars depends on our understanding of convection, 
which remains one of the key-open questions in stellar mod¬ 
elling (e.g. see Gastellani et al. 1971a; Ghiosi 2007; Salaris 
2007; Bressan et al. 2015). Crucially, stringent tests of mod¬ 
els have been limited so far by the lack of observational con¬ 
straints specific to the internal structure of evolved stars. 

In this work we propose a way forward. We argue that 
the combination of two observational constraints, i.e. the lu¬ 
minosity of the AGB bump and the RC period spacing of 
gravity modes, provides us with a decisive test to discrimi¬ 
nate between competing models of HeCB low-mass stars. 

We have computed a series of stellar models with vari¬ 
ous prescriptions for the transport of chemicals and for the 
thermal stratification of near-core regions, and using differ¬ 
ent evolution codes (MESA, BaSTI, and PARSEC). First, 
we used these models to make predictions about the dura¬ 
tion of the RC and eAGB phases, the luminosity of the AGB 
bump, and C/0 ratios in WDs. A summary of the models 
characteristics can be found in Table 1. We then focussed 
on the prediction of seismic observables (see Sec. 5). We 
found that the asymptotic period spacing of gravity modes 
depends strongly on the prescription adopted (with differ¬ 
ences up to about 40%, when comparing the BS and HOV 
model). 

We complemented this analysis by a numerical compu¬ 
tation of adiabatic oscillation frequencies. This allowed us to 
confirm that the asymptotic approximation (Eq. 2) is a good 
representation of the period spacing of gravity-dominated 
modes. Moreover (see Sec. 5.3), the detailed behaviour of 
the period spacing of g modes shows the seismic signature 
of sharp variations in the Brunt-Vaisala frequency, which 
could potentially give additional information about near¬ 
core features (localised chemical composition gradients and 
near-discontinuities in the temperature gradient). 

We then presented (Sec. 6) a first comparison between 
our predictions and the observational constraints obtained 
from the analysis of Kepler light curves (Pinsonneault et al. 
2014; Mosser et al. 2014). We found evidence for the AGB- 
bump among Aep/er targets, which allowed us to make a first 
combined analysis of classical (AGBb luminosity) and seis¬ 
mic (RC AHi) constrains. Our main conclusion is that, while 
standard models (BaSTI-SC, PARSEC-LOV) are able to re¬ 
produce the luminosity of the AGBb, they cannot describe 
satisfactorily the distribution of the observed period spacing 
of RC stars while models with high overshooting (HOV), al¬ 
though giving a much better description of the observed RC 
AHi distribution, fail to reproduce the AGBb luminosity. We 
then suggest a candidate model to describe simultaneously 
the two observed distributions: a model with a moderate 
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overshooting region in which we apply an adiabatic thermal 
stratification. This prescription (which we have tested using 
PARSEC, see PARSEC-MPC model) gives indeed a better 
description of the observations. 

At this stage of the analysis we are however pre¬ 
vented from drawing any further quantitative conclusions. 
To achieve the latter, we will follow two complementary ap¬ 
proaches. On the one hand, we will couple our models with 
TRILEGAL, to generate synthetic stellar populations which 
can be quantitatively compared with the observed compos¬ 
ite disk population. We will also investigate in detail pos¬ 
sible observational biases, both in the target selection and 
in the detection of the period spacing from oscillations spec¬ 
tra. On the other hand, to limit/quantify such biases, we will 
test our models considering stars in the clusters NGC6791 
and NGC6819, in which oscillations were detected in ~ 30 
HeCB giants (e.g., see Stello et al. 2011; Miglio et al. 2012; 
Corsaro et al. 2012). 
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Figure Al. Carbon and Oxygen profile of 1.5-M0 models at the 
first AGB-TP. 


APPENDIX A: CARBON AND OXYGEN 
CONTENT 

As mentioned in Sec. 2, an additional consequence of adopt¬ 
ing different mixing schemes in the HeCB phase is the re¬ 
sulting chemical profile of the C-O core, which affects the 
chemical profile of a WD at the beginning of the cooling 
sequence. 

To provide the reader with a complete set of predic¬ 
tions, we show in Figure Al the chemical-composition pro¬ 
file of AGB models at the first thermal pulse. The results are 
consistent with what Straniero et al. (2003) found in models 
of a 3 Mq star. They showed that the maximum of Aq for 
a bare-Schwarzschild and induced overshooting (plus semi¬ 
convection) models lies outside the former convective core 
while in models with high overshooting it corresponds to 
the value found in the former mixed core. Moreover, in our 
HPC model the Xq is higher than the HOV model, mostly 
due to the larger duration of the he-burning core phase (fig¬ 
ure Al, left panel) in which *^C(a,7)*®0 reaction has more 
time to convert carbon into oxygen. 

This paper has been typeset from a l^X/IATl^X file prepared by 
the author. 
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